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^ ■ Abstract 



KA , Experimental results on jet physics at LEP and SLC are reviewed and compared 

with perturbative QCD predictions. The discussion includes determinations of 
the strong coupling 05, measurements of event shape distributions and jet cross 
sections, studies of subjet multiplicities and tests of QCD coherence. 



* Invited talk at the 18th Johns Hopkins Workshop on Current Problems in Particle Theory, Theory 
meets experiments, Florence, September 1994. 



1. Introduction 

The large electron-positron colliders LEP and SLC have provided many detailed studies 
of hadronic final states. These include jet cross sections, direct-photon production, hadron 
spectra and, more recently, studies of fragmentation functions and their scaling violations. 
Most of these topics are discussed in several review papers (see, for instance, Refs. fll-ll). 
In this contribution, rather than presenting a general review on hadronic physics at LEP 
and SLC, I concentrate only on jet physics. The reason is the following. 

On the experimental side, LEP and SLC are remarkable machines for studying hadronic 
interactions. They operate at a very high centre-of-mass energy, a/s ~ 91 GeV, and provide 
a very large event statistics (due to the large cross section at the Z° resonance) and a very 
clean environment (as in any e~^e~ collider, the final state is not contaminated by the debris 
of initial-state hadrons). 

On the theoretical side, if we consider jet physics (i.e. high-energy hadronic observables 
which are infrared and coUinear safe) we are dealing with a very predictive theory, namely 
perturbative Quantum Chromodyanamics (QCD). We have at our disposal a well defined 
theoretical framework (with essentially only one free parameter, the strong coupling as) 
and quite accurate calculations. 

Combining these experimental and theoretical facilities, we are thus in a position to 
perform 'precision tests' of QCD. The aim of the present contribution is to show that at 
LEP/SLC we can investigate strong-interaction physics to an accuracy better than 10%, 
which is quite an achievement in this field. 

The outline of this review is as follows. In Sect. 2, 1 recall the general theoretical features 
of jet observables and their treatment within perturbative QCD. Section 3 is devoted to 
the measurements of as from fully inclusive quantities, namely the hadronic widths of the 
Z^ boson and the r lepton. More detailed tests of QCD from event shape distributions 
are considered in Sect. 4. Here, I present results based on analyses in two-loop order and 
including the all-order resummation of logarithmically enhanced contributions. A summary 
of as determinations at LEP and SLC is also included in this Section. Section 5 deals with 
the definition of jets and studies of jet cross sections. A discussion on subjet multiplicities, 
QCD coherence and properties of quark and gluon jets is presented in Sect. 6. Finally, in 
Sect. 7, I briefiy summarize the requirements for further progress in the field. 

2. Jet observables in perturbative QCD 

Jet quantities are, by definition, hadronic observables which turn out to be infrared and 
collinear safe . In other words, their actual value does not vary if the final state changes 
by the addition of one more particle which is soft or collinear to another particle. 

The main theoretical feature of jet observables is that they can be computed in QCD 
perturbation theory. This is not a trivial statement. In general, evaluating QCD Feynman 
diagrams in terms of partons (quarks and gluons), one finds integrals which are divergent 
in the low-momentum and/or small-angle regions. In the case of jet quantities, the coher- 
ent sum over different soft and collinear partonic states leads to the cancellation of these 
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divergences. As a result, jets observables are finite (calculable) at the partonic level order 
by order in perturbation theory. 

Denoting by i? a generic (dimensionless) jet variable, we can compute its corresponding 
perturbative expansion in the form 

R = l + as{Q)Ri + al{Q)R2 + ... . (1) 

The first non-trivial term (asRi) in Eq. (|lD, the second term {a'gR2), etc. represent re- 
spectively the leading-order (LO) contribution to R, the next-to-leading order (NLO) con- 
tribution and so on. The typical energy scale, at which the observable R is considered, is 
denoted by Q. 

The expansion parameter in Eq. (P is as{Q), the QCD running coupling. It is the sole 
(apart from quark masses) free parameter in the theory. QCD does not predicts its actual 
value but only its energy behaviour. Indeed, due to the property of asymptotic freedom, the 
running coupling as{Q) decreases as the energy Q increases, according to the approximate 
logarithmic behaviour as{Q) ~ [PolnQ^/A"^)'^, A being the fundamental scale of QCD. 

Asymptotic freedom is due to the non-abelian gauge structure of the theory. Thus, the 
experimental observation of the running of as{Q) according to QCD is a fundamental test 
of the underlying gauge theory. For this reason, in the past few years much effort has been 
devoted to measurements of as{Q) (see Refs. [^|6|-p!I| and Sees. 3 and 4). 

The QCD running of as implies that the effective coupling as{Q) is small at high energy 
Q. This property justifies the use of perturbation theory for predicting jet observables, at 
least at asymptotic energies. However, just because of its perturbative nature, the QCD 
running can be hidden in higher-order corrections by the replacement as{Q) = a^ [1 + 
K{Q) as{Q) + ...], Ois being the value of as at a fixed (and arbitrary) energy scale. It 
follows that a LO calculation gives only the order of magnitude of a certain observable. The 
accuracy of the perturbative QCD expansion is instead controlled by the size of the higher- 
order corrections. Any definite QCD prediction thus requires (at least) a NLO definition 
of as and a NLO perturbative calculation. 

In the following I shall consider the NLO definition of as as given in the MS renor- 
malization scheme. In this scheme the relation between as{Q) and the QCD scale Aj^ is 

m 
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where 127ij3o = 33 — 2Nf, 247r^/3i = 153 — 19Nf and Nf is the number of quarks whose 
mass is smaller than Q. 

As for the power series expansion in as beyond the LO, one has to be more precise than 
in Eq. (|^). In general, the perturbative series has the following structure 



R{as{fi),Q/l2;y) = I + as{iJ,)Ri{y) + a^iid) R2iy) - Ri{y) (3o\nQ /jj, 



+ 0{aUf^) In' Q'/f,\aUfi)) , (3) 

where Q is the typical energy scale of the process and y denotes any other kinematic scale 
involved in the definition of the observable R. 



In Eq. @ we have made explicit the dependence on the renormahzation scale /i which 
has to be introduced in any higher-order calculation in order to control (regularize) ultravi- 
olet divergences. Obviously, the observable i? is a renormahzation group invariant quantity 
if computed to ah orders in as- The /i-dependence is an artifact of the truncation of the 
perturbative expansion at a fixed order in as- This implies that a) the renormahzation 
scale dependence is formally an higher-order effect (i.e. dR/d\nfi = C(a|) in Eq. (||)) 
and, correspondingly, b) higher-order terms have an explicit ln((5//i)-dependence (see the 
0(a^g In Q'^/fx^) term in Eq. (^). However, the /i-dependence can be numerically large if i) 
the perturbative functions Ri{y), R2{y), ■- are large, and/or ii) fi is very different from the 
typical energy scale Q of the process (in this case the lnQ//i-terms in Eq. (^ are large and 
may spoil the convergence of the perturbative expansion). For these reasons, assuming a 
well-behaved perturbative expansion, fi should be set equal to Q and ^ variations (typically 
by a factor of four) around this value can be used for estimating the theoretical uncertainty 
due to uncalculated higher-order contributions. 

Before discussing in detail jet measurements at LEP/SLC, I have to add a last general 
comment about non-perturbative effects. Any hadronic obsvervable has necessarely a non- 
trivial non-perturbative component because of the hadronization of quarks and gluons. 
Therefore, in general we should write 

R = Rp,n{as{fi),Q/fi;y) + 0{{l/QY) , (4) 

where the last term represents the non-perturbative contribution. One of the main prop- 
erties of the jet observables, ensuing from their infrared and collinear safety, is that the 
non-perturbative component is suppressed by some inverse power of the energy as the en- 
ergy increases (i.e. p > 1 in Eq. (^). Therefore, in principle, the hadronization effects 
are negligible at asymptotic energies. In practice, since we are dealing with experiments at 
finite energy, we should try to estimate the size of the non-perturbative corrections. 

At present there are essentially two methods for estimating non-perturbative contri- 
butions. The first method is based on the Operator Product Expansion (OPE) [l^ and 
consists in relating the non-perturbative term in Eq. (^) to the matrix elements of some 
local operators. These matrix elements are not computable in perturbation theory but are 
universal, in the sense that their value does not depend on the particular jet quantity con- 
sidered. Therefore one can estimate them from low-energy measurements (or from lattice 
calculations) and evaluate the corresponding correction to the jet observables. 

The second method is based on Monte Carlo event generators |ll3| - |l^j . The Monte Carlo 



generators used at LEP/SLC consist of two different components. In the perturbative com- 
ponent, starting from the hard energy scale of the process, one generates a parton shower 
describing quark and gluon production according to (approximate) QCD matrix elements. 
In the non-perturbative component, at a certain scale of the order of 1 GeV, partons from 
the shower are converted into hadrons according to some phenomenological hadronization 
model. The hadronization parameters are tuned in order to reproduce the experimental 
data. Therefore, comparing jet quantities as obtained by Monte Carlo simulations both 
from partons at the end of the QCD shower and from particles after hadronization, one is 
able to estimate the size of the non-perturbative corrections. 



3. Fully inclusive quantities 

Fully inclusive quantities are jet observables depending on a single momentum scale. 
Their perturbative expansion in Eq. (^) is thus a simple power series in as with constant 
(y-independent) coefficients Rn- Moreover, due to the highly inclusive nature of these 
quantities, they are expected to be marginally affected by non-perturbative corrections 
(p > 2 in Eq. (|)). For these reasons, these jet observables are particularly suitable for a^ 
determinations. Two of them have been investigated at LEP so far: the hadronic width of 
the Z° boson and the hadronic branching ratio of the r lepton. 

3. 1 Hadronic branching ratio of the Z^ boson 

The hadronic branching fraction Rz of the Z" boson is defined by the ratio of the 
hadronic and leptonic widhts: 

p rhad(Mz) 

and is computable in terms of electroweak and QCD corrections as follows 

riep(Mz) = rr. + ri , (6) 

rhad(Mz) = Fg^ RqcD + Tew RqCD ■ iV 



Here F^ include all the known electroweak corrections |T^ within the standard model. 

The QCD contributions Rqcd ^^^ evaluated starting from the imaginary part of n^*^(Q^), 
the correlation functions of the vector [i = V) and axial {i = A) currents. These quantities 
are completely known[] up to next-to-next-to-leading order (NNLO) in perturbation theory 
(i.e. to relative accuracy 0(a|) with respect to the lowest-order approximation) |T^,pO[1, 



including all the relevant corrections due to finite quark masses [21,22 



In spite of these very accurate theoretical results, last year at the EPS Conference 
I pointed out P] that there were some discrepancies in the actual numerical implementations 
of these calculations. The effect of these discrepancies turned out to be of the same size as 
the expected theoretical uncertainties. 



This problem has been considered in a recent analysis by the authors of Ref . fl^ . They 



have shown that the numerical discrepancies are eliminated after full updating of all the 
programs |2^ . Moreover, they have provided a simple effective (factorized) formula for Rz 
and have carefully estimated its theoretical accuracy. 

The QCD corrections Rqcd ^o ^^e vector and axial part of the electroweak current are 
different. Therefore, in general, starting from Eqs. (^ and (|^) one cannot write down a 
simple power series in as for the ratio Rz in Eq. (|^). Nonetheless, in the proper range of 
electroweak parameters, one can derive an approximate factorized expression for Rz. The 

^The ©(ttg) -singlet part of the axial current, which was previously missing, has been computed recently 



effective formula obtained in Ref. [0 i^ 



Rz^R^ 



1 + 1.060— + 0.90 ( — ) -15 



^+o.9or^V-i5r^^ 



Ro = 19.943 , (8) 

wfiere as stands for as{Mz). In tfie range 0.10 < a^ < 0.15, Eq. (||) agrees witli tlie full 
perturbative formula for Rz to an accuracy better than Aas < 0.0001. 

Using the expression (j^), one can traslate theoretical uncertainties on the value of Rz 
into corresponding uncertainties on as{Mz)- The estimate in Ref. Ijl^ gives 



Aas{Mz) = ±0.002 (ew) ± 0.002 (QCD) t^m i^t, Mh) • (9) 

The electroweak error (ew) is mainly due to electroweak corrections to the Z^ -^ bb vertex. 
The QCD error is dominated by unknown higher orders in perturbation theory and is 
estimated by varying the renormalization scale n within the range Mz/4: < fi < Mz- 
Note that non-perturbative corrections to Rz are negligible because of the large value 
Mz ^91 GeV of the Z^ mass: from the OPE analysis the non-perturbative contributions 
turn out to be proportional to (l/M^)^ (even including some additional corrections of the 
type AV^I, one obtains Aas < 0.001). 

The dominant source of uncertainty in Eq. (|^) is due to the values of the masses of the 
top quark (Mt) and Higgs boson (Mh). The coefficients in Eq. (^ refer to Mt = 150 GeV 
and Mh = 300 GeV, whilst the corresponding errors in Eq. (|^) are obtained by considering 
variations in the range 200 GeV > M^ > 100 GeV, 1 TeV > M^ > 60 GeV. 

Including the data collected in the 1993 run, the updated LEP average for Rz is Rz = 
20.795 ± 0.040 P3. From this value and using Eq. (M) one gets: 



a5(Mz) = 0.124 ±0.006(exp.)+°:°°^(th.) , (10) 

where the experimental error is dominated by the event statistics. 

This result is perfectly consistent with the value as{Mz) = 0.126 ± 0.005(exp.) ± 
0.002(Mi:^) (and Mi(GeV) = 173 ± 13(exp.) ± 19{Mh)) obtained from a standard-model fit 
to the Z° lineshape and asymmetries at LEP [|2^. However the latter is more dependent 
on the parameters of the electroweak theory (in particular, Mt). 

3.2 Hadronic branching ratio of the r lepton 

An independent determination of as in NNLO can be obtained from the hadronic width 
of the T lepton. This quantity is indeed theoretically related to the current correlation 
functions 11*^*^ (Q^) via a simple momentum sum rule ||2^. The corresponding hadronic 
branching ratio is given by 

Rr = ^^I'T ""^ ^]'''^-^ = R^'^ [1 + 5pert. + 5non-pert.] , (H) 

B(r -^ Vrl^i) 



■''This formula is valid for massless leptons. Corrections due to finite lepton masses can easily be taken 
into account WA. 



where R^^^ is the parton model value and the perturbative QCD correction (as evaluated 
from the three-loop calculation of 11*^*^ (Q^), described in the previous subsection) is 

,^„,. ^ + 5.2 (^)^ 26.4 (^)'. (12) 



Unlike the case of Rz, the non-perturbative contributions are potentially sizeable at a 
scale as low as the r mass Mr = 1.78 GeV. Using the OPE and estimating the relevant 
matrix elements by means of QCD sum rules, the non-perturbative term in Eq. (0) is 
found to be [|25[ 

'^non-pert. = -0.02 ± 0.01 , (13) 



which is quite a small contribution (although of the same order as the 0(«|) term in 
Eq. (|12D, since ^^(M^) ~ 0.3). 



Although this estimate is not completely unambiguous, at present we are more confident 
on the value of (^non-pert.- As a matter of fact, the estimated value in Eq. ([T3|) has been 
found consistent with the hadronic mass spectrum of the r lepton as measured by the 
ALEPH Collaboration i6i. 



The average value^ of the results published by the four LEP experiments ||2^ is Rr 



3.617±0.034. Using this value and the theoretical predictions in Eqs. (0)-([Tj) one obtains 



asiMr) = 0.36 ± 0.02(exp.) ± 0.04(th.) . (14) 

The way of estimating the theoretical uncertainties on the value of as from r decay is 
still a matter of discussion [[f|,p8| . The main points regard the validity of the OPE close to 
the resonance region E^ and possible corrections to the (NNLO) perturbative running of 



as for scales as low as Mr [0 . The theoretical error quoted in Eq. (|T^) takes into account 



the uncertainty on (5non-pert. in Eq- (0)) the variation of the renormalization scale between 
1 GeV and 3 GeV and the effect of adding to Eq. (^) a contribution of order ±100(as/7r)^. 

As expected from QCD, the value of as{Mr) in Eq. (jT^) is significantly larger than the 
value of as{Mz) obtained from Rz- Extrapolating the result (|1^ from the r mass to the 
Z^ mass using the perturbative QCD running (see Eq. (|])), one finds 

asiMz) = 0.122 ± 0.002(exp.) ± 0.004(th.) , (15) 

where the relative size of the errors is decreased because of the logarithmic dependence of 
0:5 on the energy. 

4. Event shape distributions 

The most detailed QCD tests performed so far at e+e~ colliders are based on studies of 
shape variables and jet cross sections. 



§This average value does not include two new measurements of Rt by the ALEPH and CLEO Collab- 
orations B. These results are still preliminary and differ each other by almost three standard deviations. 



The shape variables are global jet observables characterizing the structure of the hadronic 
final states. One of these variables is the thrust T, which is defined as 

T = Max %^^ , (16) 

and thus maximizes the total longitudinal momentum (along the unit vector n) of the final- 
state particles pi in a given event. For a two-jet event we have T = 1, whilst a spherical 



event has T = 1/2. Many other shape variables can be defined [^], with the only constraint 
of being infrared and collinear safe observables. 

In the following I denote by y a generic shape variable whose two-jet limit corresponds 
to the region y —>■ (for instance, y = 1 — T). The corresponding shape- variable event 
fraction is defined by 

ry , ,1 da 
I dy--— 
/o cr ay' 

where Q is the e~^e~ centre-of-mass energy and fi is the renormalization scale. 



R{asif^),Q/f^;y)= j dy'-^j^ , (17) 



The event fraction in Eq. (0) is computable in QCD perturbation theory in the form 
of Eq. (H). Note that now the coefficients Ri{y), R2{y), ■ ■ ■ are not c-numbers but functions 
of the actual value y of the shape variable. Therefore, unlike the case of the fully inclusive 
quantities, by studying shape variable distributions one can not only measure as but also 
perfom more effective tests of the QCD matrix elements. The price which has to be paid 
is that, at comparable energy, the non-perturbative effects are stronger than for the com- 
pletely inclusive observables. Any time that the degree of inclusiveness is degraded, the 
effective power p, which controls the energy behaviour of the non-perturbative component 
(see Eq. (^), decreases. According to the naive expectation, the hadronization corrections 
to the event shapes are nominally of the order A/Q. This expectation is in agreement with 
the non-perturbative effects estimated by Monte Carlo generators and has been confirmed 
by calculations of renormalon contributions [^] . 

4-1 QCD studies to complete C(«|) 

The conventional approach for studying shape variables is based on the comparison 
between data and NLO QCD calculations. For all the relevant shape variables the 0{as) 
and C(a|) coefficients Ri{y) and R2{y) in Eq. (^ have been evaluated numerically^] by 



Kunszt and Nason |^, using the two- loop matrix elements computed in Ref. |3^. In these 
analyses the renormalization scale fi is usually set equal to Q. Note that, due to the explicit 
y-dependence of the coefficients Rniy), by varying /i one can modify the shape in y of the 
QCD predictions. Therefore, a more empirical approach is often considered: /x is left as a 
free parameter and fitted to the data together with as- 



A summary |11| of as determinations from NLO calculations at LEP and SLC |36-38 



is presented[] in Fig. 1. All these measurements are pretty consistent and the average value 

^In the case of the energy-energy correlation function, the numerical disagreement among different 
calculations pi^,pql has been solved recently [B4| in favour of Ref. [pl| . 

"The most recent (and complete) results from the SLD Collaboration are not reported in Fig. 1. How- 
ever, their inclusion does not change the average value in Eq. 



(taking into account experimental and theoretical correlations) for as is 

as{Mz) = 0.119 ± 0.001(exp.) ± 0.006(th.) . (18) 

The error is dominated by theoretical uncertainties (common to all the measurements) 
due to hadronization corrections and higher-order contributions. The hadronization effects 
are estimated via Monte Carlo event generators by comparing the corresponding results 
at parton and hadron level. The overall size of these corrections is typically between 
5 and 15% (see, for instance Fig. 2). The relative effect among different Monte Carlo 
generators is instead of the order of few percent and is usually assigned as the corresponding 
hadronization uncertainty. 

Higher-order contributions are estimated from renormalization scale variations. As a 
matter of fact, the peculiar feature of these QCD tests and as determinations to 0{a^) 
is the following |3^j3^. Far away from the two-jet region, the NLO expression (|^) with a 
renormalization scale yU ~ Q gives good fits to the data. On the contrary, reasonable fits 
in the two-jet region can be achieved only by using renormalization scales fi much smaller 
than Q (as small as few GeV!) and values of as smaller than those obtained in the multi-jet 
region (Fig. 3). 



Fig. 1: Compilation of measurements of as{Mz) from event shapes, jet rates, energy 
correlations and scaling violations, in (9(a|), at LEP and SLC [|11[]. 



Fig. 2: DELPHI data on the heavy jet mass Mh (the mass of the heavier hemisphere with 
respect to the thrust axis) compared with the 0(q;|) QCD predictions (with a renormal- 
ization scale factor / = /x^/Q^ = 0.25). The hadronization corrections apphed to the data 
are shown below the distribution. 



Fig. 3: Dependence of as{Mz) (sohd curves) and x^/d.o.f. (dashed curves) on x^ = n/Q 
for 0{a'g) fits to the OPAL data on thrust (T), heavy jet mass (Mh), total (Bt) and wide 
(Bw) jet broadening. 

The reason for this strong scale dependence of the shape- variable event fractions is that 



the corresponding perturbative functions Rn{y) in Eq. (|^) are large in the two-jet region 



(y — > 0). For instance, in the case of thrust the actual calculation gives ||39 
RTiy = l-Tl^^l-CF^\n'il-T) + Ucp^)\n\l~T) + Oial\n%l-T)). (19) 
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The double logarithmic contributions a^ ln^(l — T), a| ln^(l — T) ... are due to the brems- 
strahlung spectrum of soft and coUinear gluons. Although infrared and collinear singulari- 
ties cancel in jet observables upon adding real and virtual contributions, in the two-jet limit 
real emission is strongly inhibited. The ensuing mismatch of real and virtual corrections 
generates logarithmically-enhanced terms which spoil the convergence of the perturbative 
expansion in as- 

Observables dominated by two-jet configurations are thus affected by a large and sys- 
tematic theoretical uncertainty due to higher-order logarithmic corrections. Reliable pre- 
dictions can be obtained only computing these corrections and, if possible, resumming them 
to all orders in as- 

4-2 QCD studies using resummed calculations 

A detailed understanding of logarithmically enhanced terms now exists for many shape 
variables |40| - |4^ , namely those for which the small-y logarithms L = Inl/y exponentiate 
^ . The shape variable event fraction R{as,y) can thus be written as follows (in order to 



simplify the formulae I set /i = Q) 

R{as; y) = C{as)n<^s, L) + D{as] y) , (20) 

where 

oo 

C{as) = l + ^a«^ 

n=l 

oo n+1 

lnE(«5,l/) = EEG'n„^«5^"^ 
n=l m=l 

= Lgi{asL)+g2iasL) + asg3{asL)^ , (21) 

and D{as', y) vanishes as y — > order by order in perturbation theory. The word exponen- 
tiation refers to the fact that the terms a^L"^ with m > n+1 are absent from In R{as] y), 
whereas they do appear in R{as',y) itself. In the expression ( pOD the singular Iny depen- 
dence is entirely included in the effective form factor S. The function gi resums all the 
Zeadmg' contributions agL""*"^, while g2 contains the nea;t-to-/ea(im(jf logarithmic terms OgL^, 
and g^ etc. give the remaining subdominant logarithmic corrections a^L™ with < m < n. 



Equation (|2T|) represents an improved perturbative expansion in the two-jet region. 
Once the functions gi have been computed, one has a systematic perturbative treatment of 
the shape distribution throughout the region of y in which asL^l, which is much larger 
than the domain asL"^ -C 1 in which the as perturbative expansion (H) is applicable. 
Furthermore, the resummed expression (pT]) can be consistently matched with fixed-order 



calculations. In particular, one can consider the next-to-leading logarithmic approximation 
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(NLLA) as given by the functions gi and (72 and combine them with the 0{a'g) results 
in eq. (^ (after subtracting the resummed logarithmic terms in order to avoid double 
counting), to obtain a prediction (NLLA+0(a|)) which is everywhere at least as good as 
the fixed-order result, and much better as y becomes small. 



Fig. 4: (a) Measured distributions of thrust (r = 1 — T), heavy jet mass (p = Mfj/Q"^), 
total (-Bt) and wide (-Biy) jet broadening, compared with fits to the O^a"^) and to the 
resummed NLLA+(9(q;|) calculations, with a renormalization scale factor fi/Q = 1 in 
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both cases. The (b) hadronization and (c) detector corrections apphed to the data are also 
shown. 

Extensive experimental studies based on NLLA+(9(a|) calculations have been carried 
out during the last two years [^,|U,|^. As expected |^| from the improved theoretical 



accuracy of these predictions, it has been shown that the resummed calculations have 
a reduced dependence on the renormalization scale and, in particular, remove the need 
to choose 'unphysical' (very small) renormalization scales in the two-jet region (see, for 
instance, Fig. 4). 



Fig. 5: Compilation of measurements of as{Mz) from LEP and SLC, using resummed 
NLLA+C(a|) calculations Q. 



A partial summary of as{Mz) from resummed calculations is given in Fig. 5 [|TT|]. 
This summary does not include a very recent and detailed analysis performed by the SLD 
Collaboration [^. Combining the new SLD result {as{Mz) = 0.118 ± 0.006) with those 



of the LEP experiments, I obtain the (correlated) average value 

as{Mz) = 0.122 ± 0.002(exp.) ± 0.005(th.) . (22) 

This result is in good agreement with that obtained from analyses in 0(a|) alone 
and has a comparable uncertainty. Note however that the theoretical errors are treated 
differently. All the central values in Fig. 5 refer to the same renormalization scale value 
fj' = Q and scale uncertainties are evaluated by varying fi by (approximately) a factor of four 
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around Q. This range includes the best-fit values for /i: renormalization scales very different 
from Q are not only theoretically disfavoured but they also fail in describing the data. Non- 
perturbative effects are again estimated from Monte Carlo event generators. However, since 
modern Monte Carlo simulations operate by generating parton configurations of arbitrarily 
large multiplicity (typically much larger than the maximum of four partons involved at order 
(9(a|)), they are better suited to estimating the hadronization corrections for resummed 
calculations. In summary, using resummed calculations there is much less freedom to define 
the central value of as and its theoretical uncertainty and a more consistent QCD picture 
emerges. 

Ji^.S Summary of measurements of as at LEP and SLC 

The measurements of as at LEP and SLC can be summarized as follows (Fig. 6). We 
have two largely independent determinations of as at the scale Mz from Rz and event 
shape distributions in NLLA+0(a|). They respectively give (see Eqs. (|10D and (p2|)): 

Q;g(Mz) =0.124 ±0.008 , (23) 

as{Mz) =0.122 ±0.006 . (24) 

A third determination of as at the scale Mj. is provided by i?^: 

as{Mr) =0.36 ±0.05 . (25) 



Fig. 6: Summary of measurements of as from Rz and Rr at LEP, and from event shape 
distributions at LEP and SLC. 

The two values oias{Mz) in Eqs. (^) and (p^ are very consistent with each other, thus 
proving that perturbative QCD is extremely successful in describing hadronic interactions 
at the energy ^/s ^ 91 GeV of the Z^ boson mass. Moreover, as shown in Fig. 6, combining 
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these two determinations of as{Mz) with that of asiM^-) from r decay, we have a strong 
evidence for the QCD running of as from LEP (and SLC) data alone! Assuming the QCD 
running of Eq. (|[), the value of a5(M^) in Eq. (^) corresponds to as{Mz) = 0.122±0.005. 
Therefore, considering that the corresponding energy scales differ by almost two orders of 
magnitude {Mz — 50 Mr), the agreement among these determinations of as is a very 
significant test of QCD. 

Let me also recall that the measurements of as at LEP/SLC are in good agreement with 
those obtained from other processes. The world summary of as presented by S. Bethke 



at the Montpellier Conference [11] is reported in Fig. 7. The data points denoted by 
circles at the extremes of the plot correspond to determinations performed at LEP/SLC. 
Further evidence for the QCD running is due to the measurements at intermediate values 
of the energy Q, which are obtained mainly from deep inelastic lepton-hadron scattering 
(DIS). The values of as{Mz) from DIS are slightly smaller than those from LEP/SLC 
and lead to the consistent (and dominated by theoretical uncertainties) world average 
a5(Mz) = 0.117 ±0.006 [i|H. 



Fig. 7: Summary of measurements of as from different processes and comparison with 



perturbative QCD expectations [p]] . 

5. Jet cross sections 

A jet is qualitatively defined as a collimated spray of high-energy hadrons. It is likely to 
be produced by hard scattering of partons and thus it can be regarded as a universal signal 
of parton dynamics at short distances. However, for the purpose of performing accurate 
quantitative studies, one needs a precise definition of jet. Essentially, one has to specify 
how low-energy particles are assigned to jets, in order to have infrared finite cross sections. 
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5.1 Jet algorithms 



The standard jet definition in e^e~ annihilation p8| , pT| amounts to introducing a di- 
mensionless resolution variable yij = dij/Q"^ for every pair (i, j) of particles (jets). Then 
the particles (jets) with the minimum yij are merged until a fixed resolution ?/cut is reached 
{ijij > ycut)- The final-state particles in each event are therefore classified in a well defined 
number of jets, depending on the resolution ?/cut- The main feature of this definition is 
that the corresponding iterative procedure of clustering type provides an unambiguous and 
exhaustive assignment of particles to jets. 

Several jet clustering algorithms for e~^e~ annihilation are available [^,^. Different 
jet algorithms are specified by the definition of the dimensionful resolution variable dij. As 
a result of many thoretical and phenomenological investigations carried out in the last few 
years P3| , PD[ - ^ , the theoretically favoured resolution variable turns out to be: 



4^^ = min 2 {E^ ,E]){1- cos % ) . (26) 

Here Ei and Ej are the particle (jet) energies in the e~^e~ centre-of-mass frame and 6ij 
is their relative angle. This resolution variable reduces to the minimal relative transverse 
momentum k\^j in the limiting case of small relative angles (6'y — > 0). For this reason the 
algorithm is known as fc^-algorithmQ. 

In order to motivate the preference for the fcx-algorithm, let me comparep] its features 
with those of an older jet definition as given by the JADE algorithm |^. In the latter the 



jet resolution variable is essentially the invariant mass d\j = 2EiE j (1 — cos 9ij) ~ (Pi+Pj^ 
of the pair of particles (jets). 

The resolution variables c?L^ and d^- both vanish in the soft limit Ei,Ej -^ (the 
two algorithms are thus infrared safe) but they behave differently for small (and finite) 
particle energies. It follows that they treat soft radiation in a different way. In particular, 
since dlj depends on the product EiEj, the JADE algorithm prefers to merge soft particles 
first, even if they are far apart in angle. In the case of the /c^-algorithm, the resolution 
variable d^j-^ is instead diagonal with respect to particle energies (i.e. the product EiEj 
is replaced by Ef or Ej). Hence, soft particles are merged with the energetic particle 
closest in angle. The different jet classification is clearly shown by the L3 event reported 
in Fig. 8 IQ. One can say that the JADE algorithm induces strong attractive kinematic 
(due to the jet definition and independent of the underlying dynamics) correlations among 
soft particles. On the contrary, the fcx-algorithm, where these correlations are absent, 
leads to the several advantages: i) it avoids a non-intuitive classification of events and an 
unnatural assignment of particles to jets (soft and wide-angle jets); ii) it reduces the size of 
the non-perturbative corrections due to the hadronization process (since soft particles are 
merged with the energetic particle closest in angle, soft fragmentation products are likely 
to be assigned to the 'parent' jet); Hi) theoretical calculations in perturbative QCD are 
more reliable, because multi-parton kinematics does not dominate multi-parton dynamics 
in higher perturbative orders. 

**This algorithm was first discussed at the Durham Workshop on Jet studies at LEP and HERA, De- 
cember 1990 pi], and is sometimes referred to as the Durham algorithm. 
^^A more detailed discussion can be found in Ref. IBOll. 
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The last point is particularly evident in the small-|/cut region, where the QCD pertur- 
bative expansion of jet cross sections is dominated by large double logarithmic corrections 
of the type a§ln™j/cut (^ < Sn). These contributions, whose origin is due to the emission 
of soft and collinear gluons (as in the case of event shapes in the two-jet region), can be 
resummed to all orders in as if jets are defined using the fcx-algorithm. On the contrary, in 
the case of the JADE algorithm, the attractive kinematic correlations among soft particles 
prevent the implementation of the resummation procedure [^ . 



Fig. 8: A 3-jet event as 'seen' by (left) the JADE algorithm and (right) the fc^-algorithm. 

One more attractive feature of the fc^-algorithm is its applicability to processes in- 
volving initial-state hadrons (i.e. deep inelastic lepton-hadron scattering, photoproduction 
processes and hadron-hadron collisions). The generalization to hadron collisions of e^e~ 
clustering algorithms is by no means trivial, because one has to face the problem of deal- 
ing with the soft remnants of the incoming hadrons and factorize them from high-p^ jets 



produced by hard scattering of partons. As proposed in Ref. |5^, resolution variables 
of transverse-momentum type are particularly suitable for this purpose and, actually, a 
/c^-clustering algorithm for hadron collisions has been set up in Refs. [Q,^. Its use for 
studying jet physics at HERA and at the Tevatron collider may offer some advantages 
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54 -57 1 with respect to standard cone algorithms 5^ . 



Fig. 9: Comparison between ALEPH data on the differential two-jet rate -D2(ycut = Us) 
and resummed+(9(Q;|) calculations corrected (solid curve) for hadronization effects. 

5.2 Jet rates and multiplicities 

In the last few years, the fcx-algorithm has replaced the original JADE algorithm in 
most of the jet analyses carried out in e~^e~ annihilation, including determinations of as 
and studies of jet topology and coherence effects. 

The basic jet measurements one can consider are the n-jet rates Rn-jet = o'n-jet/o'TOT, 
defined by the ratio between the cross section an-jet for producing n jets (at the resolution 
scale i/cut) and the total cross section axoT- In Fig. 9, the ALEPH data |^ on the dif- 
ferential two-jet rate D2{ycnt) = dR2- jetijj cut) I dy cut are compared with the corresponding 
QCD predictions in 0{a'g) and including the all-order resummation of the logarithmic con- 
tributions a§ln^"?/cut and d^\v?^~^ y^nt |0- The value of as extracted from this analysis 



contributes an entry in the list reported in Fig. 5. 

As a further example of jet studies using the /cx-algorithm, let me consider the aver- 
age jet multiplicity {rijet) = Y^n>2^ Rn-jet- This measurement can be considered as an 
alternative QCD test with respect to the mean multiphcity of hadrons produced in e~^e~ 
annihilation. In the case of the hadron multiplicity one is interested in its dependence on 
the e~^e~ centre-of-mass energy (see, for instance. Sect. 6.2 in Ref. [|I|]). Here, at a fixed 
centre-of-mass energy, one can instead study the jet multiplicity as a function of the jet 
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resolution ^cut- 



{ujet) - 2 ~ Cpa 

I/cut <K 1 



TT^CAa + ^C\a^ + 



(27) 



Fig. 10: L3 data on the average jet multiplicity (njei) compared with NLLA+(9(a|) 
predictions. 

An appealing feature of {rijet) is that, unlike the hadron multiplicity, it is an infrared 
safe observable and thus its absolute normalization (and not only its dependence on the 
energy) is computable in perturbation theory. Moreover, the increase of the jet multiplicity 
is highly sensitive to non-abelian effects, i.e. to multiple jet production by gluon cascades. 
This effect is particularly evident in the small-?/cut region where the QCD calculation, 
to leading accuracy in the double logarithmic expansion parameter a = (as/^n) In ?/cut, 
predicts: 

— C^a H 

12 360 

Note that in this approximation the deviation of Eq. ( P^ ) from the double logarithmic 
behaviour as In ycut is entirely due to non-abelian (i.e. proportional to powers of the gluon 
charge Ca = Nc) contributions. The QCD predictions for (rijet) in NLLA+C(a|) |Q have 
been succesfuUy tested by the L3 and OPAL Collaborations |^ (see Fig. 10). 

5.3 Other jet studies 

The gauge structure of QCD is completely determined by the fact that the gauge group is 
SU{Nc) with Nc = 3 colours, and that the quarks belong to the fundamental representation. 
This fixes completely the relative magnitude of all the couplings of the theory. Roughly 
speaking, to the splitting processes q —^ qg, g -^ gg and g — > qq, one can associate splitting 
probabilities respectively proportional to Cp = {N^ — l)/2Nc, Ca = Nc and Tr = 1/2. 
Events with more than four final-state jets in e'^e" annihilation give access to all the 
relevant splitting processes and thus multijet events offer the possibility of performing very 
detailed tests of the gauge structure of the theory. 

Angular correlations within four-jet events have been investigated by the LEP Col- 
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laborations |5^. In these investigations, the values of the colour factors Cp, Ca, T^ in 
the two-loop matrix elements [^] are left as free parameters and fitted to the LEP data. 
The results of the fits are consistent with the colour gauge group S'f/(3) whilst the only 
other groups with three quark colours (namely, 5*0 (3) and the abelian toy model intro- 
duced in Ref. [Q) are excluded. Unfortunately, since only C(a|) matrix elements are at 



present available for these four-jet studies, the ensuing calculations are just LO predictions. 
Therefore it is very difficult to assess the theoretical accuracy of these analyses. 

Many more studies of jet properties have been performed at LEP and SLC Among 
them, the results, which have been recently reported [^,^| on the evidence for differences 



between quark and gluon jets, are particularly relevant. This and other jet topics are 
reviewed elsewhere ||l|,|^,§]. In the next section, I discuss in some detail a study of the 
internal structure of jets based on the measurement of the subjet multiplicity inside jets. I 
think this measurement is quite interesting for addressing issues like QCD coherence and 
properties of quark and gluon jets. 

6. Sub-jet multiplicity and QCD coherence 

6. 1 Quark and gluon jets 

One of the main features of perturbative QCD is the different colour charge carried by 
quarks and gluons. The bremsstrahlung spectra for the emission of a soft gluon from a 
hard quark and from a hard gluon differ only in their relative normalization. The emission 
probability from a quark (or an antiquark) is proportional to the square |Tgp = Cp = 4/3 
of the (non-abelian) quark chargeQTg, whilst the emission probability from a gluon is pro- 
portional to the square |Tgp = Ca = 3 of the gluon charge Tg. Among the quantities 
which are easily measurable in the experiments, the average particle (parton) multiplicity 
is the most sensitive to the radiation of soft gluons. If we compare the particle multiplic- 
ities Mg and Afq in a gluon and a quark jet, we may thus expect a ratio Afg/Afq which is 
approximately unity at low energies and increases with the energy. At asymptotic energies 
(in practice, as soon as Mg,J\fq ^ 1) this ratio should be equal to the ratio of the (squares 
of the) colour charges, Afg/Afq ~ Cf/Ca = 9/4, i.e. larger than a factor of two. 

This naive expectation has two main problems, besides possibly large subasymptotic 
corrections |^3|]. One problem is related to the hadronization effects. The hadron multi- 



plicity is not an infrared safe observable and thus it cannot be computed in QCD pertur- 
bation theory. In other words, we have poor theoretical control over the non-perturbative 
contribution to the ratio of the hadron multiplicities: some systematic difference in the 
hadronization of quark and gluon jets could significantly contribute to a deviation from the 



asymptotic multiplicity ratio |^ . 



The other problem is related to the definition and identification of quark and gluon jets. 
As discussed at the beginning of Sect. 5, the quantitative definition of jet is a highly non- 
trivial issue, which can be solved only by introducing jet-finding algorithms. The problem 
is amplified in the case of a distinct definition of quark and gluon jets because, although we 



*The quark and gluon colour charges T, and Tg are SU{?>) colour matrices respectively in the funda- 
mental and adjoint representations. 
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observe jet-like objects, quarks and gluons are not experimentally accessible. Thus, quark 
jets and gluon jets can be singled out by comparing different processes in which only quarks 
or gluons contribute at the level of the naive parton model. Alternatively, in a given process 
one can anti-tag a gluon jet by identifying heavy-quark decays [^,^]. No matter how quark 
and gluon jets are defined, their definition is theoretically biased by our underlying partonic 
picture. More quantitatively, we can say that quark and gluon jets cannot be defined in 
a universal way to an accuracy better than 0{as), whether due to the dependence on the 
process or the dependence on the observable. Moreover, although nominally of 0{as), this 
dependence can be numerically large because typically enhanced by logarithmic coefficients 
(as discussed in Sect. 5.1, even the double logarithmic contributions a^ln^ycut to the jet 
rates do depend on the jet definition). 



Some years ago |^J], we propose a method aimed at investigating differences between 
quark and gluon jets by overcoming the problems discussed so far. The idea is to select two- 
jet and three-jet event samples with the /i;_|_-algorithm at a fixed value of the jet resolution 
parameter ?/cut = Vi- Then the fc^-algorithm is again applied by clustering subjets with a 
smaller resolution j/cut = 2/o < 2/i and counting the number of subjets in each jet sample. 

Within this approach the hadronization corrections are under better control because 
hadrons are replaced by subjets, which are infrared and collinear safe objects. Moreover, 
e"'"e~ annihilation is a point-like source of qq events which undergo fragmentation through 
the splitting process qq -^ qqg. Therefore, comparing the two-jet and three-jet samples, 
one has access respectively to the properties of a mixture of qq jets and qqg jets without 
further specification and identification of quark and gluon jets. 

In particular, one can consider the ratio of the average subjet multiplicities M3 and M2 
in three-jet and two-jet events. By fixing yi and decreasing i/q, the ratio M3/M2 is expected 
to vary in the range: 

3 ^ M3 ^ Af, + Af^ + ^g ^ 2Cf + Ca ^ 17 

2- M2^ Uq+Uq ~ 2Cf 8 ■ ^ ' 

The value 3/2 on the left-hand side of Eq. (|^) corresponds to yo = yi and is simply due to 
kinematics (at yo = yi subjets and jets coincide and there are exactly 3 jets in the three- 
jet sample and 2 jets in the two-jet sample). The right-hand side of Eq. (^) is instead 
approached at |/o -C ?/i- Actually, the value 17/8 is that naively expected as yo/yi ^- on 
the basis of the (asymptotic) charge counting rule Afg/Afg = Afg/Afg ^ Ca/Cf- 



The theoretical [^ and experimental [^] analyses show that this naive expectation 
is not only numerically violated but completely misleading from a physical viewpoint. In 
any physical process, jets are not produced independently because of colour conservation. 
This leads to coherence effects that, in the case of the subjet multiplicities, are responsible 
for a value of the ratio M3/M2 which remains well below 3/2 (which is the value at the 
kinematic boundary and not the asymptotic value!) for most of the yo range 



The perturbative QCD predictions obtained in Ref. |^| include the resummation of 
the leading and next-to-leading logarithmic contributions a^L^" and a^L"^"^'^ {L standing 
for both Li = —Inyi and Lq = —Inyo) to all orders in as- The actual calculation is 
quite involved, so that, in the following, I limit myself to presenting the final results and 
to providing their qualitative interpretation. 
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6.2 Two-jet and three-jet events: Arj-rapidity profiles 

Considering jets produced at the e+e~ centre-of-mass energy Q, let me introduce the 
transverse-momentum scales Qi = ^Jy\Q and Qo = y/VoQ- Because of the definition (^) 
for the resolution variable, jets selected by the fc^-algorithm at ?/cut = Hi have a natural 
angular size Oj ~ 2Qi/ Ej, Ej being the jet energy (Fig. 11a). This angular size identifies 
the jet core, i.e. the fragmentation region of the jet: all the subjets (or particles), which 
are inside a cone of aperture 9j around the jet axis, belong to the jet. However, each jet J 
contains also all the subjets which are relatively soft (with energy much smaller than Ej) 
and are produced at an angle 9 (with respect to the jet axis) in the region 9j^29^ 9jj', 9jj' 
being the angular distance to the nearest jet J' (Fig. 11). This region can be called 'interjet' 
region since its angular extension depends on the topology of the multijet configuration. 

Let me thus consider the (pseudo-)rapidity profile dN/drj of a jet, that is, the number 
of subjets as a function of their pseudorapidity rj = — lntan(6'/2) with respect to the jet 
axis. 

A two-jet event with resolution ?/cut = Vi consists of the subjets radiated by a leading 
quark and a leading antiquark (Fig. 11a). In the fragmentation region of the quark jet the 
rapidity profile increases up to a saturation value and then produces a rapidity plateau 
in the interjet region (Fig. 12). The situation is perfectly symmetric in the antiquark 



hemisphere and the result for the subjet multiplicity in the two-jet event is |Q: 



M2iQo,Qi;Q) = 2\Af,{Qo,Qi) + CpH{Qo,Q,)\n^\ . (29) 




The first term on the right-hand side of Eq. (^) represents the 'intrinsic' multiplicity of 
the quark and antiquark jets, coming from their natural angular regions. This contribution 
is indeed obtained by integrating the rapidity profile in Fig. 12 over the fragmentation 
regions. 

The subjet multiplicity ^fq{Qo,Ql) is equal to I {njetiVcnt = Ql/Ql)), {njetiVcut = 
Ql/Q"^)) being the total jet multiphcity discussed in Sect. 5.2. As a matter of fact, if 
yi = 1 (i.e. Qi = Q — Ej/2) the fragmentation region extends up to r^ = and all the 
e^e~ events are identified as two-jet events. 

The second term on the right-hand side of Eq. (^) gives the extra interjet mul- 
tiplicity between the jet cores. Note that each jet contributes to a rapidity plateau 
whose length and height are respectively ?7i = \yi2Ej/Qi ~ InQ/Qi and Cp H{Qo,Qi) ~ 
dN'q{Qo,Qi)/d\nQi. The height of the plateau is proportional to the square |Tgp = Cp 
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of the colour charge of the leading quark. 



Fig. 11: (a) A two-jet event in e^e~ annihilation. The subjets ki and k2 are produced 
respectively in the fragmentation region (I) and in the interjet region (II). (b) Typical 
kinematic configuration of a three-jet event: (A), (B), and (C) denote the angular regions 
described in the text. 



Fig. 12: Rapidity profile of a quark (antiquark) jet in a two-jet event: (I) fragmentation 
and (II) interjet regions (r^o = \n.2Ej/QQ ~ InQ/Qo, Vi = ln2_Ej/(5i ^ InQ/Qi). 

Let me then consider three-jet events. When t/i ^ 1 (as appropriate for a calculation 
to logarithmic accuracy in Inyi) a three-jet event has the typical kinematic configuration 
in Fig. lib. An antiquark (quark) jet recoils from a quark (antiquark) and a gluon jet 
which are produced at a small relative angle 6qg in the opposite hemisphere. The rapidity 
profile of the antiquark jet is similar to that in the two-jet sample. The situation is instead 
different in the hemisphere with two jets. One has to examine subjet production in three 
different angular regions (Fig. lib): (A) at small angle 6 with respect to the quark-jet axis 
{9 < 9qg); (B) at small angle 9 with respect to the gluon-jet axis {9 < 9gg); (C) at large 
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angle with respect to both the quark-jet and the gluon-jet axis {6 > 9, 



ig) 



Fig. 13: Rapidity profiles for the angular regions (A), (B) and (C) (see text) of the 
quark+gluon jet hemisphere in three-jet events {rjqg = — lntan6'qg/2). 

The angular regions (A) and (B) include the jet cores and lead to rapidity profiles 
(Fig. 13) which are similar to that observed in two-jet events. The only difference between 
the quark jet (A) and the gluon jet (B) is that the latter has a larger multiplicity Mg{Qo, Qi) 
coming from its fragmentation region]] and, correspondingly, a higher rapidity plateau. The 
height of the rapidity plateau of the gluon jet is CaH{Qo,Qi) ~ dAfg{Qo,Qi)/d\nQi, i.e. 
proportional to the square jT^p = Ca of the gluon charge. 

Subjets produced in the interjet region (C) also contribute to a rapidity plateau (Fig. 13). 
However, since their emission angle 6 is much larger than 6qg, they cannot resolve two dis- 

TThe gluon multiplicity J^g{Qo,Q) is equal to one half of the total jet multiplicity one could measure 
in quarkonium decay into two gluons. 
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tinct jets. The height of the plateau is thus proportional to the square |Tg + Tgp of the total 
colour charge of the quark and gluon jets. Because of colour conservation Tg + Tg = —Tq 
{Tq being the colour charge of the antiquark jet recoiling in the opposite hemisphere), and 
hence, the height of the plateau is proportional to |Tq + T^p = | — T^p = Cp ! 



Fig. 14: Ar^-rapidity profiles for (a) two-jet and (b) three-jet events (Lj = —Inyi, i = 0, 1). 
The suppression at small A?] in (b) is due to QCD coherence. 

The results of this discussion on three-jet events can be combined in the rapidity profile 
of Fig. 14b. Here At] = T]—\n{2Ej/QQ), where r] is the subjet pseudorapidity with respect to 
the axis of the parent jet and Ej is the energy of the parent jet. Therefore the A?7-rapidity 
profile is experimentally measurable without identifying quark and gluon jets. Integrating 
the Ay^-rapidity profile in Fig. 14b, one obtains [|6^] : 



M3(Qo, Qi; Q) = M2(Qo, Qi; Q) + Ar,(Qo, Qi) + Ca H{Qo, Qi) (r/)^ . (30) 

The term M2 on the right-hand side of Eq. (^) comes from the contribution of the 
antiquark jet and from the contributions of the angular regions (A) and (C) (which add 
themselves to reconstruct an entire quark-jet rapidity profile as in Fig. 12) in the oppo- 
site hemisphere. The remaining terms on the right-hand side of Eq. (|30|) are due to the 
additional gluon jet. Like the quark and antiquark jets, the gluon jet provides an intrin- 
sic and an interjet contribution. However, the latter is restricted to the angular region 
9qg > 9 > Qi/Ej around the gluon jet axis. Thus the length of the pseudorapidity plateau 
is Inq/Qi, where q = m.m{Eq9qg, Eg9qg) is the transverse momentum at which the gluon 
and quark jets are resolved. In Eq. (|30|) , this length is averaged over q {{rj) = (Ing/Qi)) 
with the appropriate three-jet cross section. 

In Fig. 14a, I have also reported the Ar^-rapidity profile of the two-jet sample, in order 
to make easier the comparison with the corresponding profile of the three-jet sample. We 
can see that the additional gluon jet contribution in Fig. 14b, although asymptotically 
higher by the canonical factor Ca/Cf = 9/4, is always shorter than those of the quark jets 
if Hi <^ 1. In this case the angle between the quark and the gluon jet is very small and the 
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gluon jet is strongly squeezed. This suppression of the A?7-rapidity profile in the central 
region is due to QCD coherence, that is, to destructive interference of the quark and gluon 
jets. In this large-angle region, the quark and gluon jets act coherently with a colour-charge 
factor of |Tg + Tgp = I - T^p = Cp, which is much smaller than |Tgp + \Tg\^ = Cf + Ca, 
the colour-charge factor that one would consider in the case of independent fragmentation 
of the two jets. 

A general lesson one can learn from this discussion is that quark and gluon jets do 
not live independently. Because of QCD coherence (colour conservation) their properties 
depend on their production mechanism and on the actual selection procedure. 

6. 3 Multiplicities in two- and three-jet events 



Figure 15 shows results on the ratio M3/M2 from LEP ||B^. As a consequence of the 
squeezing of the gluon jet due to QCD coherence, this ratio is predicted to fall off near the 
kinematic limit yo — yi [0. This prediction is in good agreement with the data. Also 
shown is the prediction of a toy model which does not take into account coherent effects. In 
this toy model, where jets contribute according to their energies and colour charges without 
interference, the multiplicity ratio rises monotonically towards its asymptotic value, in 
disagreement with the experiments. As a further check on the origin of the decrease in the 
ratio M3/M2, one can measure its slope at yo = yi as a function of yi. According to QCD 
coherence, by increasing yi the angle 6qg increases, thus reducing the region in which the 
quark-gluon destructive interference takes place. The slope of M3/M2 is thus predicted to 
decrease at large yi, in agreement with the data ||65[ . 



Fig. 15: L3 data on the ratio M3/M2 compared with the QCD predictions in Eqs. (|2 
(solid curve) and with the predictions of an incoherent toy model (dotted curve). 

The measurement of the subjet multiplicity in two-jet and three-jet events can be def- 
initely considered a detailed test of (perturbative) QCD coherence. Since it is based on 
infrared and collinear safe observables which do not require any identification of the gluon 
jet, this test of coherence is certainly more 'model independent' than the celebrated string 
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effect |66l-|68i. 



Fig. 16: OPAL data on the subjet multiplicity of two-jet events compared with different 
Monte Carlo models. The fractional difference between models and data is also shown. 



The perturbative QCD predictions for M3/M2 are confirmed by the LEP data up to 
yo — 5 • 10^^. This value of yo corresponds to a minimal transverse momentum Qq — 2 GeV 
at LEP energies. It is thus evident that perturbative QCD fails at smaller values of yo 
because the hadronization effects dominate. On the contrary, Monte Carlo event generators 
can be used for investigating the development of a jet by varying y^ from the formation scale 
(1/0 = Vi) to the scale at which individual hadrons are resolved as subjets {y^ ~ 10^^ at Q = 
Mz). The measurements of subjet properties performed at LEP have been compared with 
different Monte Carlo models |6^. These include generators like JETSET 

HI 



131, HERWIG 



ARIADNE [|l5f 



which implement coherence effects in the parton shower evolution, 
and generators like CO JETS |jT6|, which are based on independent fragmentation. As 
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shown by the representative study in Fig. 16, all the models describe well the data. This 
is not surprising: the various (hadronization) parameters have been tuned to describing 
many hadronic observables and such a tuning can mask the dynamics differences among 
the models. However, it is interesting to note that, although the agreement with the data 
is always very good for yo ^10~^ (i.e. at scales where hadronization effects dominate), the 
Monte Carlo generators with no QCD coherence perform worse at the true subjet level 
{yo^lO~^). Typically, the incoherent models produce too much radiation and too many 
subjets in this phase space region. This analysis is a further evidence in favour of coherence 
and shows that the study of subjets can be a useful tool for investigating the transition 
from the perturbative to the non-perturbative phase of QCD. 

7. Outlook 

As discussed in the previous sections, with the advent of the high-energy and high- 
statistics experiments at LEP/SLC and of accurate theoretical calculations, 'high-precision' 
tests of QCD have become possible. These tests include measurements of as and its energy 
dependence, detailed analyses of two-jet and three-jet events and studies of QCD coherence. 

Further improvements in the field require mainly progress in the theory. Future 0{a^) 
calculations for event shapes and jet cross sections can be very useful to check the consis- 
tency of the present NLO analyses [^|. Accurate studies of multijet (n > 4) final states can 



be performed only if the corresponding QCD predictions in NLO become available. The 
increasing capabilities of tagging 6-quark jets by micro-vertex detectors allow to carry out 
precise experimental studies in this field, but complete NLO calculations are still missing 

This list does not exhaust the demands to the theory. One of the main problems indeed 
remains our poor understanding of the hadronization process and, in general, of the non- 
perturbative region. In order to deal with this problem we need better theoretical methods 
and more phenomenological investigations. In this respect, studies of infrared and collinear 
safe observables at the boundary of the perturbative region can be of some help. As briefly 
discussed at the end of Sect. 6, the small-yo behaviour of the subjet multiplicity in two- 
jet and three-jet events is an example of this kind of studies. Further calculations and 
predictions for observables which are suitable for similar investigations are warranted. 

Acknowledgements 

I wish to thank Yu.L. Dokshitzer, F. Fiorani, M. Olsson, L. Trentadue, G. Turnock and 
B.R. Webber for a long collaboration on jet physics. I am also grateful to S. Banerjee, 
S. Bethke, P.N. Burrows, G. Cowan, T. Hebbeker, R. Miquel, M. Schmelling, D.R. Ward 
and A. Wehr for many useful discussions. This research is supported in part by the EC 
Programme "Human Capital and Mobility" , Network "Physics at High Energy Colliders" , 
contract CHRX-CT93-0357 (DC 12 COMA). 

References 



27 



1. T. Hebbeker, Phys. Rep. 217 (1992) 69. 

2. S. Bethke, in Proc. of the Aachen Conf. QCD - 20 Years Later, eds. P.M. Zerwas and 
H.A. Kastrup (World Scientific, Singapore, 1993), pag. 43; preprint PITHA 94/29. 

3. J. Drees, in Proc. of the Aachen Conf. QCD - 20 Years Later, eds. P.M. Zerwas and 
H.A. Kastrup (World Scientific, Singapore, 1993), pag. 106. 

4. B.R. Webber, plenary talk at the 27th Int. Conf. on High Energy Physics, Glasgow, 
July 1994 (quoted as ICHEP 94 in the following), preprint Cavendish-HEP-94/15. 

5. G. Sterman and S. Weinberg, Phys. Rev. Lett. 39 (1977) 1436. 

6. S. Bethke and S. Catani, in Perturbative QCD and Hadronic Interactions, Proc. 27th 
Rencontres de Moriond, ed. J. Tran Thanh Van (Editions Frontieres, Gif-sur-Yvette, 
1992), pag. 203. 

7. G. Altarelli, in Proc. of the Aachen Conf. QCD - 20 Years Later, eds. P.M. Zerwas 
and H.A. Kastrup (World Scientific, Singapore, 1993), pag. 172. 

8. S. Bethke, in Proc. of the 26th Int. Conf. on High Energy Physics, Dallas, August 
1992, ed. J.R. Sanford (AIP, New York, 1993), pag. 81. 

9. S. Catani, in Proc. of the Int. Europhysics Conf. on High Energy Physics, HEP 
93, Marseille, July 1993, eds. J. Carr and M. Perrottet (Editions Frontieres, Gif-sur- 
Yvette, 1994), pag. 771. 

10. Particle Data Group, Review of Particle Properties, Phys. Rev. D 50 (1994) 1297. 

11. S. Bethke, talk at the Conference QCD 94, Montpellier, July 1994, preprint PITHA 
94/30. 

12. K. Wilson, Phys. Rev. 179 (1969) 1499; K. Symanzik, Comm. Math. Phys. 32 
(1971) 49; C. Callan, Phys. Rev. D 5 (191972) 3302; R. Brandt, Fortschr. Phys. 
18 (1970) 249; M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B147 
(1979) 385, 448, 519. 

13. T. Sjostrand, Computer Phys. Comm. 39 (1986) 347, Computer Phys. Comm. 43 
(1987) 367; M. Bengtsson and T. Sjostrand, Nucl. Phys. B289 (1987) 810. 

14. B.R. Webber, Nucl. Phys. B238 (1984) 492; G. Marchesini and B.R. Webber, Nucl. 
Phys. B310 (1988) 461; G. Marchesini et al.. Computer Phys. Comm. 67 (1992) 465. 

15. U. Pettersson, preprint LU TP 88-5; L. Lonnblad and U. Pettersson, preprint LU TP 
88-15; L. Lonnblad, preprint LU TP 89-10. 

16. R. Odorico, Computer Phys. Comm. 32 (1984) 139, Computer Phys. Comm. 59 (1990) 
527; P. Mazzanti and R. Odorico, Nucl. Phys. B370 (1992) 23. 

17. T. Hebbeker, M. Martinez, G. Passarino and G. Quast, Phys. Lett. 331B (1994) 165. 

18. K.G. Chetyrkin and O.V. Tarasov, Phys. Lett. 327B (1994) 114; S.A. Larin, T. van 
Ritbergen and J. A.M. Vermaseren, Phys. Lett. 320B (1994) 159. 

28 



19. S.G. Gorishny, A.L. Kataev and S.A. Larin, Phys. Lett. 259B (1991) 144; L.R. Sur- 
guladze and M.A. Samuel, Phys. Rev. Lett. 66 (1991) 560. 

20. B.A. Kniel and J.H. Kiihn, Phys. Lett. 224B (1989) 229, NucL Phys. B329 (1990) 

557. 

21. K.G. Chetyrkin and A. Kwiatkowski, Phys. Lett. 305B (1993) 285; K.G. Chetyrkin, 
Phys. Lett. 307B (1993) 169; K.G. Chetyrkin and J.H. Kiihn, Phys. Lett. 308B (1993) 
127; K.G. Chetyrkin and A. Kwiatkowski, Phys. Lett. 319B (1993) 307. 

22. K.G. Chetyrkin and J.H. Kiihn, preprint TTP 94-08. 

23. G. Burgers, W. Holhk and M. Martinez, program BHM; W. Holhk, Fortschr. Phys. 
38 (1980) 3, 165; D. Bardin et al., program ZFITTER, version 4.6, preprint CERN- 
TH 6443/92; G. Montagna et al., program TOPAZO, Phys. Commun. 76 (1993) 
328. 

24. D. Schaile, plenary talk at ICHEP 94. 

25. E. Braaten, S. Narison and A. Pich, Nucl. Phys. B373 (1992) 581; F. Le Diberder 
and A. Pich, Phys. Lett. 289B (1992) 165. 

26. ALEPH Coll., D. Buskulic et al, Phys. Lett. 307B (1993) 209. 

27. E. Fernandez, in Proc. of the Int. Europhysics Conf. on High Energy Physics, 
HEP 93, Marseille, July 1993, eds. J. Carr and M. Perrottet (Editions Frontieres, 
Gif-sur-Yvette, 1994), pag. 705. 

28. S. Narison, preprint CERN-TH.7188/94; A. Pich, talk at the Conference QCD 94, 
Montpellier, July 1994. 

29. T.N. Truong, preprint EP-CPTh.A266-1093 (EP-CPTh.A266-1093-REV). 

30. G. Aharelli, preprint CERN-TH. 7246/94. 

31. Z. Kunszt, P. Nason, G. Marchesini and B.R. Webber, in 'Z Physics at LEP 1', CERN 
89-08, vol. 1, p. 373. 

32. B.R. Webber, preprint Cavendish-HEP-94/7. 

33. A. Ah and F. Barreiro, Phys. Lett. 118B (1982) 155, Nucl. Phys. B236 (1984) 269; 
D.G. Richards, W.J. Stirhng and S.D. Ellis, Phys. Lett. 119B (1982) 193, Nucl. Phys. 
B229 (1983) 317; N.K. Falk and G. Kramer, Zeit. Phys. C42 (1989) 459. 

34. E.W.N. Glover and M.R. Sutton, preprint DTP/94/80. 

35. R.K. Ellis, D.A. Ross and A.E. Terrano, Nucl. Phys. B178 (1981) 421. 

36. L3 Coll., B. Adeva et al., Phys. Lett. 248B (1990) 464, Phys. Lett. 257B (1991) 469; 
ALEPH Coll., D. Decamp et al, Phys. Lett. 255B (1991) 623, Phys. Lett. 257B (1991) 
479; DELPHI Coll., P. Abreu et al., Zeit. Phys. C54 (1992) 55; OPAL Coll., P.D. 
Acton et al., Zeit. Phys. C55 (1992) 1; S. Bethke and J.E. Pilcher, Ann. Rev. Nucl. 
Part. Sci. 42 (1992) 251. 

29 



37. Mark-II Coll., S. Komamiya et al., Phys. Rev. Lett. 64 (1990) 987; DELPHI Coll., P. 
Abreu et al., Phys. Lett. 311B (1993) 408; OPAL Coll., R. Akers et al., Zeit. Phys. 
C63 (1994) 197. 

38. SLD Coll., K. Abe et al., Phys. Rev. Lett. 71 (1993) 2528, preprints SLAC-Pub-6451 
(1994), SLAC-Pub-6641 (1994). 

39. P. Binetruy, Phys. Lett. 91B (1980) 245. 

40. S. Catani, G. Turnock, B.R. Webber and L. Trentadue, Phys. Lett. 263B (1991) 491. 

41. S. Catani, G. Turnock and B.R. Webber, Phys. Lett. 272B (1991) 368, Phys. Lett. 
295B (1992) 269. 

42. J. Kodaira and L. Trentadue, Phys. Lett. 123B (1982) 335; J.C. Collins and D.E. 
Soper, Nucl. Phys. B197 (1982) 446, Nucl. Phys. B284 (1987) 253; G. Turnock, 
preprint Cavendish-HEP-92/3; R. Fiore, A. Quartarolo and L. Trentadue, Phys. 
Lett. 294B (1992) 431. 

43. S. Catani, Yu.L. Dokshitzer, M. Olsson, G. Turnock and B.R. Webber, Phys. Lett. 
269B (1991) 432. 

44. S. Catani, Yu.L. Dokshitzer, F. Fiorani and B.R. Webber, Nucl. Phys. B377 (1992) 
445. 

45. S. Catani, G. Turnock, B.R. Webber and L. Trentadue, Nucl. Phys. B407 (1993) 3. 

46. ALEPH Coll., D. Decamp et al., Phys. Lett. 284B (1992) 163; L3 Coll., O. Adriani 
et al., Phys. Lett. 284B (1992) 471; OPAL Coll., P.D. Acton et al., Zeit. Phys. C59 
(1993) 1; DELPHI Coll., P. Abreu et al., Zeit. Phys. C59 (1993) 21. 

47. L3 Coll., S. Banerjee, in Proc. of the Int. Europhysics Conf. on High Energy Physics, 
HEP 93, Marseille, July 1993, eds. J. Carr and M. Perrottet (Editions Frontieres, 
Gif-sur-Yvette, 1994), pag. 299; TOPAZ Coll., Y. Ohnishi et al., Phys. Lett. 313B 
(1993) 475; TPC/Two-Gamma Coll., D.A. Bauer et al., preprint SLAC-Pub-6518 
(1994); TOPAZ Coll., M. Aoki et al., contributed paper GLS0944 to ICHEP 94. 

48. JADE Coll., W. Bartel et al., Zeit. Phys. C33 (1986) 23. 

49. S. Bethke, Z. Kunszt, D.E. Soper and W.J. Stirhng, Nucl. Phys. B370 (1992) 310; N. 
Brown and W.J. Stirling, Zeit. Phys. C53 (1992) 629. 

50. S. Catani, in QCD at 200 TeV, Proc. 17th Eloisatron Project Workshop, eds. L. 
Cifarelli and Yu.L. Dokshitzer (Plenum Press, New York, 1992), pag. 21. 

51. Yu.L. Dokshitzer and M. Olsson, Nucl. Phys. B396 (1993) 137. 

52. Report of the Hard QCD Working Group, in Proc. Durham Workshop on Jet Studies 
at LEP and HERA, J. Phys. G17 (1991) 1537. 

53. N. Brown and W.J. Stirling, Phys. Lett. 252B (1990) 657. 

54. S. Catani, Yu.L. Dokshitzer and B.R. Webber, Phys. Lett. 285B (1992) 291. 

30 



55. S. Catani, Yu.L. Dokshitzer, M.H. Seymour and B.R. Webber, Nucl. Phys. B406 
(1993) 187; S.D. Ellis and D.E. Soper, Phys. Rev. D 48 (1993) 3160. 

56. M.H. Seymour, in QCD and High Energy Hadronic Interactions, Proc. 28th Rencon- 
tres de Moriond, ed. J. Tran Thanh Van (Editions Frontieres, Gif-sur-Yvette, 1993), 
pag. 141; S. Catani, in QCD and High Energy Hadronic Interactions, Proc. 28th 
Rencontres de Moriond, ed. J. Tran Thanh Van (Editions Frontieres, Gif-sur-Yvette, 
1993), pag. 257; B.R. Webber, J. Phys. G19 (1993) 1567. 

57. M.H. Seymour, Zeit. Phys. C62 (1994) 127, Nucl. Phys. B421 (1994) 545. 

58. J. Huth et al., in Proc. of Research Directions for the Decade, Snowmass 1990, ed. 
E.L. Berger (World Scientific, Singapore, 1992), pag. 134. 

59. OPAL Coll., M.Z. Akrawy et al., Zeit. Phys. C49 (1990) 49; L3 Coll., B. Adeva et 
al, Phys. Lett. 248B (1990) 227; ALEPH Coll., D. Decamp et al., Phys. Lett. 284B 
(1992) 151; DELPHI Coll., P. Abreu et al., Zeit. Phys. C59 (1993) 357, contributed 
paper GLS0180 to ICHEP 94; OPAL Coll., contributed paper GLS0597 to ICHEP 
94. 

60. G. Kramer, Springer Tracts in Modern Physics (1984) 102. 

61. OPAL Coll., P.D. Acton et al., Zeit. Phys. C58 (1993) 387; DELPHI Coll., K. 
Hamacher et al., contributed paper GLS0187to ICHEP 94; ALEPH Coll., contributed 
paper GLS0539 to ICHEP 94; SLD Coll., Y. Iwasaki and U. Tohoku, SLAC-Pub-6597 
(1994). 

62. ALEPH Coll., contributed paper GLS0545 to ICHEP 94. 

63. A.H. Mueller, Nucl. Phys. B213 (1983) 85 (Erratum Nucl. Phys. B241 (1984) 141), 
Nucl. Phys. B228 (1983)) 351; Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan, in Per- 
turbative Quantum Chroniodynamics, ed. A.H. Mueller (World Scientific, Singapore, 
1989), pag. 241. 

64. S. Catani, B.R. Webber, Yu.L. Dokshitzer and F. Fiorani, Nucl. Phys. B383 (1992) 
419. 

65. L3 Coll., S. Banerjee, L3 Note 1171 contributed paper to the 26th Int. Conf. on High 
Energy Physics, Dallas, August 1992; ALEPH Coll., D. Decamp et al.. Note A845 
contributed paper to the 26th Int. Conf. on High Energy Physics, Dallas, August 
1992; L3 Coll., O. Adriani et al., Phys. Rep. 236 (1993) 1; OPAL Coll., R. Akers et 
al., preprint CERN-PPE/94-52. 

66. JADE Coll., W. Bartel et al., Phys. Lett. lOlB (1981) 129. 

67. B. Andersson, G. Gustafson and T. Sjostrand, Phys. Lett. 94B (1980) 211, Zeit. Phys. 
C6 (1980) 235; Ya.I. Azimov, Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan, Phys. 
Lett. 165B (1985) 147. 

68. DELPHI Coll., P. Aarnio et al., Phys. Lett. 240B (1990) 271; OPAL Coll., M.Z. 
Akrawy et al., Phys. Lett. 261B (1991) 334; L3 Coll., B. Adeva et al, Zeit. Phys. C55 
(1992) 39. 

31 



69. New Techniques for Calculating Higher Order QCD Corrections, Proc. ETH Work- 
shop, Zurich 1992, ed. Z. Kunszt, preprint ETH-TH/93-01. 

70. A. Ballestrero, E. Maina and S. Moretti, Phys. Lett. 294B (1992) 425. 



32 



